Short-time Fourier transform laser Doppler holography 
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We report a demonstration of laser Doppler holography at a sustained acquisition rate of 250 
Hz on a 1 Megapixel complementary metal-oxide-semiconductor (CMOS) sensor array and image 
display at 10 Hz frame rate. The holograms are optically acquired in off-axis configuration, with a 
frequency-shifted reference beam. Wide-field imaging of optical fluctuations in a 250 Hz frequency 
band is achieved by turning time-domain samplings to the dual domain via short-time temporal 
Fourier transformation. The measurement band can be positioned freely within the low radio- 
frequency spectrum by tuning the frequency of the reference beam in real-time. Video-rate image 
rendering is achieved by streamline image processing with commodity computer graphics hardware. 
This experimental scheme is validated by a non-contact vibrometry experiment. 



Though effective for single-point analysis [l|, laser 
Doppler measurements are more difficult to perform in 
wide-field imaging configuration, because of a technolog- 
ical challenge : digital image frames have to be read out 
at kHz rates and beyond to perform short-time discrete 
Fourier transforms (DFT) [2j. Recently, image-plane 
laser Doppler recordings with a high throughput CMOS 
camera in conjunction with short-time DFT calculations 
by a field programmable gate array (FPGA) reportedly 
enabled continuous monitoring of blood perfusion about 
the mm/s. Full- field flow maps of 480 x 480 pixels 
were rendered at a rate of 14 Hz, obtained from image 
recordings at a frame rate of 14.9 kHz Q. For transient 
dynamics imaging of faster phenomena, high throughput 
laser Doppler schemes were designed by multipoint [4| or 
time multiplexing [H[ approaches. High speed holography 
enabled offline vibrometry from time-resolved optical 
phase measurements [f|. Heterodyne holography, as a 
variant of time-averaged holography 0, HJ with a strobe 
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;9] or a frequency-shifted reference beam (lol 
appropriate for steady-state (at the scale of the exposure 
time) mechanical vibrations mapping. Advances in 
reconstruction techniques of optically-measured digital 
holograms with Graphics Processing Units (GPUs) 
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1 31 ] have led to real-time holographic screening of a 
sing le vibration frequency, demonstrated in this regime 



In this letter, we report an experimental demonstra- 
tion of video-rate image reconstruction and display of 
laser fluctuation spectra from high speed holographic 
measurements. Sustained Fresnel reconstruction of 
off-axis holograms at 250 Hz and 10 Hz rendering by 
short-time DFTs is performed. Images and RF spectra 
of a thin metallic plate's out-of-plane vibration modes 
around 3.2 kHz are presented. 

The optical setup is similar to the one reported in the 
demonstration of video-rate vibrometry at a single fre- 



quency 1J], at the difference that a high throughput 



CMOS camera is used to achieve megapixel recordings 
at 250 frames per second. An off-axis, frequency-shifted 
Mach-Zehnder interferometer is used to perform a mul- 
tipixcl heterodyne detection of an object field E beating 
against a separate LO field -Elo- i n reflective geometry. 
The main optical radiation field is provided by a 100 
mW, single-mode laser (wavelength A = 532 run, optical 
frequency is L = 5.6 x 10 14 Hz, Oxxius SLIM 532). The 
optical frequency of the LO beam is shifted by an arbi- 
trary quantity Av in the low RF range by two acousto- 
optic modulators (AA-electronics, MT80-A1.5-VIS). The 
object studied is a a thin metallic plate with hexagonal 
holes, shined over ~ 30 mm x 30 mm with ~ 50 mW 
of impinging light. It is excited with a piezo-electric 
actuator (PZT, Thorlabs AE0505D08F), vibrating sinu- 
soidally, driven at 10 V. The structure's vibrations pro- 
voke a local phase modulation <f> (cqfTJ) of the backscat- 
tered optical field E. Interference patterns are measured 
with a Basler A504k camera (Micron MV13 progressive 
scan CMOS sensor array of 1280 x 1024 pixels, quan- 
tum efficiency ~ 25 % at 532 nm). The camera is run 
in external trigger mode at v§ = ojs/(2tt) = 250 Hz, at 8 
bit /pixel quantization. Images of the central 1024 x 1024 
pixels region are recorded. The image acquisition is inter- 
faced with a National Instruments NI PCIe-1433 frame 
grabber. Each raw interferogram digitally acquired at 
time t, noted X(t) = \E(t) + E^o(t)\ 2 is dumped to a 
1024 x 1024 x 1 byte frame buffer in the GPU RAM of 
a NVidia GTX 580 graphics card by a CPU thread (fig. 
QJ. The object field of complex amplitude E is noted 



E = £ exp (iuht + i(f>(t)) 



(1) 



where = 27tz/l and (f>(t) is the fluctuating phase, as a 
result of optical path length modulation. The acousto- 
optic modulators enable the optical LO field of complex 
amplitude £lo to be detuned by Av = Auj/(2tt) 



Elo = £lo exp (iujht + iAwt) 



(2) 
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Holographic image rendering from each recorded inter- 
ferogram is performed with a numerical Fresnel trans- 
form. The hologram /, back-propagated to the object 
plane, is calculated by forming the Fast Fourier Trans- 
form (FFT) T of the product of I with a quadratic 
phase map, depending on the relative curvature of the 
wavefronts of E and -Elo m the sensor plane. This cal- 
culation is handled by the GPU (thread #1, Fig. [I]), 
by an algorithm elaborated with Microsoft Visual C++ 
2008 and NVIDIA's Compute Unified Device Architec- 
ture (CUDA) software development kit 3.2. The prac- 
tical implementation of free-space propagation with a 
discrete Fresnel transform 14j, [l5| yields complex- valued 



holograms / = \E\ 2 + I^loI 2 + E*E W + EEl Q recon- 
structed in the object plane. In off-axis configuration 
16], the zero-order terms \E\ 2 and £lo| 2 an d the twin- 
image term E*E^o can be filtcred-out. After filtering, 
the remaining complex- valued contribution to the off- axis 
hologram is 

H(t) = EEl Q = ££l Q exp(i<i>(t) - iAuit) (3) 

The heterodyne spectrum of the radiation field E is de- 
tected by a short-time discrete Fourier transform (DFT) 
of H(t) over N = 250 consecutive samples (fig. Q] thread 
#2), 10 times per second. The m-th Fourier component 
of the DFT, 

N 

H m {t) = S~]H(t- n/u s ) cxp (-2mmn/N) (4) 
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FIG. 1: Algorithmic layout of holographic rendering. Raw 
interferograms are recorded by the main CPU thread. Spa- 
tial Fresnel transforms are performed by a first GPU thread. 
Short-time temporal Fourier transforms are performed by a 
second GPU thread. 
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FIG. 2: Sketch of the spectral screening region. 



is a heterodyne measure of the laser fluctuation spectrum 
at time t, at frequency Av + v m 



H m (t) = H{t,Av + v m ) 



(5) 



iThe discrete frequencies v m of the measured spectra lie 
wthin the Nyquist limits of the camera bandwidth ±z^s/2, 
while the LO detuning frequency Av can be set arbitrar- 
ily by the acousto-optic modulators (fig. [5]). 

We assessed the thin metal plate's out-of-plane vibra- 
tion modes around 3.2 kHz with the presented holo- 
graphic approach. The metallic structure was excited 
sinusoidally at one (P = 1) or two (P = 2) frequencies 
, and ism 2 ■ I n cither case, the resulting out-of-planc 
motion at a given point of the surface of the plate is 



z(t) = ^z p sin (uj Mp t) 
P =i 



(6) 



where wm p = 27tzz-m p and z p are the angular frequency 
and the local amplitude of each component, respectively. 
The phase modulation of the backscattered light is 



p 

p=i 



(T) 



where <p p = Airz p /\. The temporal part of the field un- 
dergoing sinusoidal phase modulation can be decomposed 
in a basis of Bcssel functions with the Jacobi- Anger iden- 
tity. If P = 1, the object field is 



E — £ exp (ioj\J) 



E 

n— — oo 



J n (</>i)exp(mw Ml (8) 



where J„ is the Bcssel function of the first kind of rank n. 
The only remaining low frequency, camera-filtered term 
in eq. [3] beating around Ais, within the Nyquist domain 
(i.e. for \u Ml - Ai/\ < is 

H hF (t) ^ £l £Ji((j) 1 ) c^p (iuj Ml t ~ iAut) (9) 

This modulated hologram yields a single component in 
the short-time DFT spectrum H m (t), at the frequency 
v m = I'M! — Av. In the first part of the movie (media 1, 
when P = 1), the excitation frequency i^Mi was swept 
from 3210 Hz to 3290 Hz, and the LO was detuned by 
Av = 3200 Hz; the measurement frequency v m of the 
short-time DFT was swept concurrently from 10 Hz to 
90 Hz, in 5 Hz steps. The reported spectra result from 
the magnitude \H{t,Av + v m )\ averaged within the red 
square superimposed on the vibration maps. 
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EalTTI yields the frequency comb observed in figEtg), 
whose peaks are separated by \vm 2 ~ v Mi\ = 5 Hz. In 
the second part of the movie reported in media 1, the 
first excitation frequency was set to Uu 1 = 3290 Hz, 
the second one was swept from vm 2 = 3290 Hz to 
vm 2 = 3250 Hz, in 5 Hz steps. The frequency comb 
broadened with \vm 2 ~ v M t | in accordance with equations 
ITT1 and IT2"1 For a larger frequency difference \vm 2 — vmi |, 
only two lines of the comb are visible (figlSIc,f)). 




-1 25 



125 



FIG. 3: Short-time Fourier transform vibration maps and 
spectra of a metallic structure excited sinusoidally. No ex- 
citation (a). Excitation at P = 2 frequencies : um 1 — 
3200 + 100 Hz and vu 2 = 3200 - 70 Hz (b-c,e-f). Exci- 
tation at P = 2 frequencies : vm 1 = 3200 + 50 Hz and 
via 2 = 3200 + 45 Hz (d,g). LO detuning : Av = 3200 Hz; 
The movie of the experiment is reported in media 1. 



The metallic structure was then excited sinusoidally at 
two frequencies : z/Mi , and vm 2 ■ The object field under- 
going phase modulation from a double excitation (P = 2 
in eq. [7|) takes the form 

2 oo 

E = £ exp (iui-^t) J ^2 J n (4> p ) exp (inuj Mp t) (10) 

p— 1 n— — oo 

The terms of eq. [3] modulated at frequencies within the 
camera bandwidth ±ojg/2 arc actually measured. The 
others arc filtcrcd-out. The temporal part of the camera- 
filtered hologram reduces to the low frequency compo- 
nent 



H hF (t) = £l £c- lAut Cn+l,lC-„,2 + c_ nilCn+1 , 2 (ll) 
where 
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In conclusion, we performed laser Dopplcr imaging 
from sustained sampling of 1 Mega pixel interferograms 
at a throughput of 250 Mega bytes per second, and ren- 
dering of 0.25 Mega pixel off-axis heterodyne holograms 
by short-time discrete Fourier transform with a refresh- 
ment rate of 10 Hz. This demonstration was made with 
commodity computer graphics hardware. We reported 
video-rate optical monitoring of out-of-plane vibration 
amplitudes in a frequency band of 250 Hz, shifted by 3.2 
kHz from DC. This demonstration opens the way to high 
bandwidth laser Doppler holography in real time. 



Cn, P = Jni't'p) exp (inu)M p t) 
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